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INTRODUCTION 

Deep inelastic scattering (DIS) has proven to be an invaluable source of information 
about the structure of the nucleon[l], providing insight into parton distribution functions. 
Semi-inclusive deep-inelastic scattering (SIDIS) has extended our understanding further, 
allowing further access to the transverse structure of nucleons. 

With a deeper understanding of the fragmentation functions, these studies will con- 
tinue to provide us with valuable information on nucleon structure. Fragmentation func- 
tions are an important theoretical tool in the investigation of scattering reactions. This 
has led to the development and study of the NJL-jet model [2, 3, 4]. The NJL-jet model 
builds on the Field-Feynman Quark-jet model (FFQJM)[5], by using an effective chiral 
quark model framework in which calculations of both quark distribution and fragmenta- 
tion functions can be performed without introducing ad hoc parameters. Pion fragmenta- 
tion functions in the NJL-jet model were calculated in Ref. [2]. The NJL-jet model was 
then extended to include strange quark contributions, and kaon fragmentation functions 
were obtained in Ref. [3] and in Ref. [4] the fragmentation functions to vector meson 
and nucleon-antinucleons were calculated. 

FRAGMENTATION FUNCTIONS AND THE NJL-JET MODEL 

In Fig. 1, it is shown how a quark can produce a cascade of hadrons, producing what is 
interpreted as jet events in Deep-Inelastic Scattering (DIS). It is important to note that 
within the model the emitted hadrons do not interact with the other hadrons produced 
in the quark jet. An integral equation for the quark cascade process shown in Fig. 1 is 
derived in the NJL-jet model of Ref. [2]. The integral equation for the total fragmentation 
function is: 
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FIGURE 1. Quark Cascade 



The first term on the right hand side is the driving function and represents the prob- 
ability of creating a meson m carrying momentum fraction z to z + dz from the initial 
quark(i.e. hi = m). The second term represents the probability of creating the meson, m, 
further in the quark decay chain(i.e. h-2 — m). The method used to solve this equation 
approximates integrals as sums over discretized values of z, so that D'^(z) and can 
be written as vectors and the integral term (without D™(z)) can be written as a matrix, 
where the elements are their values at z = Zi- In index form the elements of the vectors 
can be written as D™- = // + gijD™j. Writing in vector form and rearranging to solve for 

D™, we obtain D™ = (I — g)~ l f, where / is taken to be the driving function and g is the 
matrix form of the integrand without the fragmentation function. 

DIHADRON FRAGMENTATION FUNCTIONS 

The extension to dihadron fragmentation functions in the NJL-jet model is now consid- 
ered. Dihadron fragmentation functions D q l! 2 {z\,Z2) correspond to the probability of 
observing two hadrons, h\ and hi, with light-cone momentum fractions z\ and Z2, re- 
spectively, fragmenting from the initial quark q. An illustration of how a quark cascade 
can produce two observed hadrons in the NJL-jet model is shown in Fig. 1, where h\ 
and h2 are the observed hadrons. 

The integral equation for the dihadron fragmentation function D q u 2 (zi,Z2) con- 
structed from Field and Feynman's Eqs (2.43a)-(2.43d) of Ref. [5] is: 

D^(zuz 2 ) = #(zO-P=^ + #fe)-^^ (2) 
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where d q l (z) and dq(rj) are the elementary splitting functions of the quark q to the 
corresponding hadron h and quark Q. The first term corresponds to the probability of 
producing hadron hi from quark q at the first step in the cascade, followed by hadron h 2 
either directly afterwards or further down the cascade chain. Similar for the second term 
with hadrons hi and h2 switched in order. The third term corresponds to the probability 
of both hadrons being produced after the first step of the cascade. 

Collecting the first two terms of Eq.(2) together, a term analogous to the driving func- 
tion of the single hadron fragmentation function is found for the dihadron fragmentation 
function. Due to the discretization method used to solve the equation, the division by T] 
causes the values of Dg' h2 (^, ^) obtained to not correspond to the discrete values in 



the chosen region for z\ and z 2 - To fix this issue, the values of z\ and Z2 were discretized 
in the region [0,1] and new variables were set as %\ = ^ and ^2 = ^. Values for t,\ and 
<^2 were discretized in the same way as z\ and zi in the region [0,1]. Rewriting the inte- 
gral part of Eq.(2) in terms of these new variables, the dihadron fragmentation function 
becomes: 
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Mathematica was used to solve for the single hadron and dihadron fragmentation 
functions. The number of points used for the dihadron cases was 200, while the single 
hadron case was 500 points. Fig. 2 shows the results for the dihadron fragmentation 
function Df n "(zi,z 2 ) solved for 50, 100, 150 and 200 points at z\ = 0.1 (left) and 
Z\ = 0.5(right). For the z\ = 0.5 case there is very good agreement for each choice 
of the number of points. For the z\ =0.1 case, sufficient convergence of the function is 
obtained by 200 points, and so this is the number of points chosen to be used. To solve 
for the delta function in the integral term the values of z\ and zi are looped over first, 
and then £1 is looped over(zi and Z2 share the same range of values as £1 and £2, and 
are discretized in the same way). This leaves ^2 with a value obtained from the delta 
function that may not be one of the discretised values. The values of the DFFs at ^2 
were obtained using linear interpolation from neighbouring discrete values. 




.« 0.02F 



0. 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



FIGURE 2. Dihadron fragmentation functions at the model scale for h\ 
number of z points at (left) z\ =0.1 and (right) z\ = 0.5 
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RESULTS 

Using the data obtained from solving for the dihadron fragmentation functions, it is 
possible to compare the solutions (z 2 ^V 2 ) to the driving function (the first two terms 
of Eq. 2). In Fig. 3 it can be seen that the driving functions contribute almost all of the 
DFF. The figure on the left shows that for an initial up quark, there is a slight difference 



between the DFF and the driving function at low values of zi for a low value of zi(=0.1). 
The strange quark has a zero value driving function, and therefore is completely made 
up of the integral term. The difference between the DFFs and the driving functions for 
each case is of the same order of magnitude as the strange quark DFF. This difference 
is only visible for the initial up quark case when the driving function becomes small, 
which occurs at low values of z\- 
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FIGURE 3. Dihadron fragmentation functions at the model scale for h\ = n + , hi = n at (left) z\ =0.1 
and (right) z\ = 0.5 



CONCLUSIONS 



In this document, results have been presented for dihadron fragmentations calculated 
within the NJL-jet model. The difference between the driving functions and the full 
DFFs was shown to become apparent only at low zi ■ The integral term for each case is 
of the same order of magnitude, though they are not equal. The difference is visible for 
the initial up quark case only where the driving function becomes small. 
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